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Using recently developed Lanczos technique we study finite-temperature properties of the 2D Kondo 
lattice model at various fillings of the conduction band. At half filling the quasiparticle gap governs 
physical properties of the chemical potential and the charge susceptibility at small temperatures. 
In the intermediate coupling regime quasiparticle gap scales approximately linearly with Kondo 
coupling as A qp /J ~ 0.3. Temperature dependence of the spin susceptibility reveals the existence 
of two different temperature scales. A spin gap in the intermediate regime leads to exponential 
drop of the spin susceptibility at low temperatures. Unusual scaling of spin susceptibility is found 
for temperatures above T c > 0.6 J. Charge susceptibility at finite doping reveals existence of heavy 
quasiparticles. A new low energy scale is found at finite doping. 
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I. INTRODUCTION 

The Kondo lattice model is one of the simplest two- 
band lattice models of correlated electrons. It is widely 
used to model heavy fermion materials where weakly in- 
teracting electrons in wide bands coexist with almost lo- 
calized electrons in unfilled orbitals of actinide or rare- 
earth elements. In heavy fermion materials a remarkable 
variety of different phases can be found at low tempera- 
tures: paramagnetic metal with large quasiparticle mass, 
anti-ferromagnetic and ferromagnetic phases, unconven- 
tional superconductivity, etc. In most of this cases, 
strong electron correlations represent the key ingredient 
of the theory that explains the rich variety of physical 
phenomena. 

In this work we investigate the Kondo lattice model, 
defined on small two dimensional square lattices with 
periodic boundary conditions. We use the finite- 
temperature Lanczos method Q . The model can be writ- 
ten as 
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<ij>s 
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where Si = J2 SS ' c i S > a s's c is an d summation < ij > runs 
over nearest neighbors. There are two distinct types of 
degrees of freedom in this model: free electrons described 
by c operators and localized spins described by S. In the 
limit when J = the two systems are decoupled which 
leads to a large degeneracy of the states due to nonin- 
teracting spins. At finite Kondo coupling J ^ the 
two systems interact. It is believed, that the interplay 
between the two degrees of freedom represents the most 
important physical mechanism of the heavy Fermion ma- 
terials. Associated with the two distinct systems are two 
competing interactions that govern the low-temperature 
physics. Finite Kondo coupling leads to formation of 
Kondo spin-singlets between the conducting electrons 
and the localized spins which screens the moments of 



localized spins. Singlet formation competes with band 
propagation of electrons in the conduction band. Cou- 
pling between kinetic energy of band electrons and lo- 
cal Kondo coupling leads to formation of the Ruderman- 
Kittel-Kasuya-Yoshida (RKKY) interaction between lo- 
calized spins. The Kondo screening on the one hand and 
the RKKY interaction on the other are in many cases 
competing interactions. Conditions, under which one or 
the other prevails depend mostly on the strength of the 
Kondo coupling, the electron filling and the dimension- 
ality of the system. 

A number of theoretical approaches has been applied 
to the investigation of the Kondo lattice model g]. In 
the strong coupling regime perturbation theory can be 
applied Large- Nf expansion can be used in the case 
of large localized spin-degeneracy Nf Q. Slave-boson 
approach Q, Gutzwiller variational treatments j5j and 
recently developed strong coupling method || have been 
successful in predicting the heavy mass of the quasiparti- 
cles, the phase diagram and the properties of the spectral 
functions. 

Numerical calculations have been mostly limited to 
one dimensional systems where various well developed 
techniques are available and finite-size effects can be 
easily controlled. Calculations on small systems have 
demonstrated that at half-filling the one dimensional 
Kondo lattice model is a spin-liquid with a finite spin 
gap . Density-matrix numerical-renormalization group 
(DMRG) calculations ||, provided accurate determi- 
nation of the spin and charge gaps as a function of 
the Kondo coupling [p|jic[|. Recently, a powerful finite- 
temperature DMRG method pl| , p^| have provided re- 
liable results for thermodynamic |13 15 and dynamic 
properties [jl5 16 1 of the model at T > 0. 

While there are many reliable numerical results of the 
Kondo lattice model in ID, much less is known about the 
model in two or three dimensions. Based on theoretical 
considerations conceptually different physical behavior is 
expected in higher dimensions. Spin-charge separation 
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exists in ID models as a consequence of strong correla- 
tions. It is reflected in different energy scales that govern 
the low-energy behavior of spin and charge excitations 
leading to a difference between the spin and the charge 
velocities. Luttinger liquid parameters define power-law 
behavior of correlation functions in ID while in higher 
dimensions exponential behavior of correlation functions 
is expected unless long-range order exists. The lack of 
long-range order in ID is responsible for Kondo screening 
to overcome the RKKY interaction for any finite J. In 
two dimensions there exists a critical value of J c /t ~ 1.4 
JL7| p"9[ below which RKKY interaction prevails and sys- 
tem orders antiferromagnetically in the case when the 
conduction band is half filled. 

The main purpose of this work is to explore thermo- 
dynamic properties of the 2D Kondo lattice model. We 
focus our investigations to intermediate and high tem- 
peratures and try to identify various energy and temper- 
ature scales that govern the spin and the charge response 
of the system. Due to small system size we are not able to 
explore extremely low temperatures since in this regime 
strong finite-size effects emerge. For the same reason 
we limit our calculations to intermediate and strong cou- 
pling regime where physics is sufficiently local so that our 
results remain valid even in the thermodynamic limit. 



II. RESULTS 

Our numerical calculations are performed by recently 
developed finite-temperature Lanczos method p| . We in- 
vestigate square lattices of N = 8 and 10 sites. Most of 
the results presented are for the N = 10 case. Stan- 
dard zero-temperature exact diagonalization results on 
small clusters are generally plagued by strong finite-size 
effects. Performing calculations at finite temperatures 
and within the grand-canonical ensemble gives us not 
only the thermodynamical properties of the system, but 
most importantly diminishes finite-size effects for T > T* 
jjj. This temperature depends primarily on the number 
of low lying excited states in the system. The T* can be 
very small when the system possesses either: a) a large 
number of low-lying energy states, or b) if physics is suf- 
ficiently local. Local physics is expected at large Kondo 
coupling where the size of the Kondo singlet is of the or- 
der of a few lattice spacings. We present results for half- 
filled conduction band case n c = 1 and at finite doping 
6 defined by n c = 1 — S. Due to particle-hole symmetry 
only 6 > is considered. In this work we restrict calcu- 
lations to thermodynamic quantities as are the chemical 
potential /i, spin susceptibility \s —< S^ a 2 > /T, charge 
susceptibility \c — —d6/dfi — {{n 2 ) — (n c ) 2 )/T and the 
specific heath c v — —Td 2 F/dT 2 . Despite small system 
size we can compute thermodynamic quantities at any 
doping S simply by choosing the appropriate chemical 
potential m. 



A. The chemical potential fj, 

In Fig. (|]a) we show chemical potential /j, as a function 
of temperature for small doping values and J/t = 2.5. 
Due to particle-hole symmetry the relation n(6 = 0) = 
is valid at any temperature. At high temperature and 
finite doping the chemical potential approaches a simple 
expression (faint dotted lines), calculated by the finite- 
temperature expansion 



M = T log 



1-6 _ _5_ f iJ 2 
1 + 6 ~ 2T 2 l~8~ 



8t 



(2) 



While the first term is model independent, the second 
term represents the first nontrivial finite-temperature 
correction. Numerical results start deviating significantly 
from the hight-temperature expression below T/t ~ 2. 
At small temperatures chemical potential approaches a 
finite value even in the limit 6^0. Infinitcsimally small 
doping away from the half-filled conduction band leads 
to an abrupt change of the chemical potential. Such be- 
havior indicates formation of the quasiparticle gap A qp 
also found in the Id DMRG calculations [Q. To inves- 
tigate the quasiparticle gap in more detail, we present in 
Fig. (|ljb) fi/J at fixed doping for different choices of the 
Kondo coupling strength on a system of N = 10 sites. We 
find that the quasiparticle gap increases almost linearly 
with the Kondo coupling J in the intermediate coupling 
range, i.e. for J/t = 1.6, 2.0 and 2.5. Extrapolated values 
of the quasiparticle gap in the limit T — ► are presented 
in table (Q) for systems of N = 8 and 10 sites. Results, 
obtained from the two systems agree reasonably well in 
the intermediate coupling range. The lack of finite-size 
effects in the limit when T — > is attributed to the exis- 
tence of the quasiparticle gap. Near the strong coupling 
limit, J/t = 10, the extrapolated values for the quasi- 
particle gap agree well with the strong coupling result 
A qp = y-2t+ ft 2 1 J = 5.72 @. 
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FIG. 1. Chemical potential (iasa function of temperature 
T; a) at different dopings S and fixed J, b) at fixed doping S 
and different couplings J/t. 



B. Spin and charge susceptibilities 

1. Strong coupling (atomic) limit 

At large values of Kondo coupling, J/t —10, physics of 
the Kondo model becomes local. We support our claim 
by results shown in Fig. (||) where we present compar- 
ison of spin (and charge) susceptibility for J/t = 10 
and analytical results obtained within the atomic limit 
of the Kondo model. In this limit the grand canoni- 
cal sum can be calculated as Z = Y^j=i e ~ <3i ~ Ej ~^ N ', 
where (3 is the inverse temperature and only 8 states 
are taken into account: the singlet state with the energy 
—3 J/4, the three- fold degenerate triplet state 
J/4 both containing one conduction electron 
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and four-fold degenerate states Es=i/2 = consisted of 
an empty and a doubly occupied conduction level each 
of them with two different spin configurations. Values of 
spin, quasiparticle and charge gap are in this limit given 
by A s /J = 1, A qp /J = 3/4, and A c /J = 1.5. At zero 
doping (8 = 0) spin and charge susceptibility are given 
by simple expressions 



Xs=P 



1 + 2 e -^/ 4 



4 + 3e-/"/4 + e 3/?J/4 ' 
4 

4 + 3e-/"/4 + e 3/3J/4 ' 



(3) 
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where (3 is the inverse temperature and fi — 0. Sus- 
ceptibility, Eq. (Q) presented in Fig. (0) follow 1/T law 
at high temperatures. At lower temperatures we see a 
peak, marked by arrows. We introduce two tempera- 
ture scales that mark the peak positions: T s /J = 0.453 
and TqpJ J — 0.386 for spin and charge susceptibility re- 
spectively. Given values were obtained analytically. At 
low temperatures both susceptibilities approach zero at 
(5 = which is consistent with the existence of a gap 
in the excitation spectrum. Low-temperature behavior 
is in both cases given by Xc.s oc (3 Exp [— 3/3J/4] which 
leads us to a conclusion, that the quasiparticle (smallest) 
gap A q p/J = 3/4 governs the low-temperature behavior 
of both susceptibilities. This is possible since a quasi- 
particle excitation modifies the charge configuration and 
also changes the spin quantum number by ±1/2. Even 
though both susceptibilities share a common gap, they 
reach a maximum at slightly different temperatures T s 
and T qp which is due to different nature of the excitation 
spectra above the gap. 

We see that numerical calculations at large J/t = 10, 
also presented in Fig. (0), agree reasonably well with an- 
alytical results at zero and finite doping. Agreement at 
finite doping is somewhat surprising since in the atomic 



limit only a single site of the Kondo lattice is taken into 
account. To understand the divergence of spin and charge 
susceptibilities at low temperatures and finite doping, we 
perform low-temperature expansion which gives us 
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Divergence at low temperatures is in the strong coupling 
regime a consequence of a degenerate level system. At 
even lower temperatures, numerical results should satu- 
rate towards finite values in both cases, the spin and the 
charge susceptibility. In the later case a deviation from 
the predicted 1/T law can be seen in Fig. (|^b) while in 
the former case such deviation is expected at even lower 
temperatures. 
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FIG. 2. Spin a) and charge b) susceptibilities \ a , \ c vs. 
T/J calculated numerically at J/t = 10 (full lines). Dashed 
curves represent results obtained in the atomic limit. Posi- 
tions of the peaks T a in a) and T qp in b) in the atomic limit 
and zero doping are indicated with arrows. 



2. Zero doping 

In Fig. (||a) we present spin susceptibilities Jxs(T/J) 
at zero doping, <5 = for different values of the Kondo 
interaction J. At high temperature (T > J) numerical 
results agree with the high-temperature expansion, 



Xs 



S-5 2 
8T 



1 - 



I - 8 2 J 



8 2 2T 



(8) 



performed to the first nontrivial order. In the intermedi- 
ate temperature regime we find rather surprising result. 
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All curves merge on a single curve for T/J > T c /J ~ 
0.6. It could be argued that this is because the high- 
temperature result in Eq. (|[) scales with J, i.e. the func- 
tion Jxs{J/T) is independent of J. However, the agree- 
ment with the Eq. ||) is only within 10% up to Tj J ~ 1 
(see Fig. (||a)) while the overlap of susceptibilities calcu- 
lated for a wide range of J/t is within a few percents. 
At low temperature, spin susceptibility reaches a maxi- 
mum at T = T s and then approaches zero. In the strong 
coupling limit spin gap is larger than the quasiparticlc 
gap however at smaller J/t this is no longer true. In the 
region of small J/t ~ 1.4 spin gap approaches zero due 
to formation of AFM order p9| . Low temperature be- 
havior of spin susceptibility in the intermediate coupling 
regime is thus governed by the spin gap. There are two 
possible approaches to estimate the spin gap using our 
method. In the zero-temperature approach the spin gap 
equals the energy difference between the lowest (S = 0) 
and the first excited (S = 1) state. At finite tempera- 
ture, the spin gap is roughly proportional to the position 
of the peak in \s (T) given by the activation temperature 
T s . We believe, that the second method, even though in- 
direct, gives results that are closer to the thermodynamic 
limit. 
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FIG. 3. Spin a) and charge b) susceptibilities Jxs, Jx? 
vs. T/J at zero doping. Faint dashed line in a) represent 
high-temperature expansion result, Eq. (^). Legends, given 
in a) apply also for b) and the inset in b). In the inset in b) 
all six different curves for J/t — 1.2 ... 10 are presented. All, 
except for J/t — 10 appear as a single line. Legend FT in a) 
indicates analytical result in Eq. (|J). 

Values of T s are presented in Table (|). As seen from 
Fig. (^a) with increasing coupling J, low-temperature 
peak in spin susceptibility moves toward higher values 
of T/J. T s therefore does not scale linearly with J for 
small J as does the quasiparticle gap (see also Table (|)). 
This points towards a non-linear dependence of spin-gap 



vs. J/t that is found in one-dimensional system [^) and 
also recent calculations in the 2D system [^9|. A qual- 
itatively good agreement is also found between T s and 
the spin gap obtained by the projection quantum Monte 
Carlo simulations G3]. 

At small values of the Kondo coupling J/t < 1.4 gap- 
less AFM long range order develops on an infinite lattice 
as a consequence of RKKY interaction pj^ , ^9[ . Uniform 
spin susceptibility in this case saturates around the tem- 
perature which is given by the RKKY interaction be- 
tween localized spins. Our results in this regime become 
less reliable at low temperature due to strong finite-size 
effects. 

The charge susceptibility Xc shows in sharp contrast 
to the strong coupling limit in many respects differ- 
ent behavior than the spin susceptibility. Starting from 
the high-temperature limit we first present the high- 
temperature result for the charge susceptibility 



Xc = 



1-5 2 
2T 



1 - 



1 



8T 2 



3J 2 



8r 



(9) 



Note that in contrast to spin susceptibility, the first non- 
trivial correction in the inverse temperature is 1/T . 
The agreement of numerical results with the high- 
temperature expansion strongly depends on J. In 
Fig. (|]b) we show Jxc(T/J) at zero doping and a wide 
range of J. At large J/t = 10, where J is the domi- 
nant energy scale in the system, numerical results agree 
with the analytical result, given by Eq. (^), down to 
T/J = 0.4. At smaller J high-temperature limit is 
reached above T/J > 2. 

The charge susceptibility is governed by a single energy 
scale, i.e. the quasiparticle gap A qp . This is reflected in 
nearly perfect scaling of J% C (T/J) for the N = 8 sys- 
tem in the intermediate coupling region 1.2 < J/t < 2.5 
(see the inset of Fig. (||b)). Scaling is due to the fact 
that the quasiparticle gap scales nearly linearly with J, 
i.e. A qp ~ 0.3 J in this regime. The quasiparticle gap 
remains finite even at small J ~ 1.2 where the spin gap 
disappears. The scaling does not persist up to the strong 
coupling limit, J/t = 10 due to a crossover regime where 
the spin gap becomes larger than the quasiparticle gap. 
The location of the peak at T — T qp seen in the tem- 
perature dependence of the charge susceptibility curves 
matches the value of the quasiparticle gap T qp = A qp ob- 
tained from measurements of the doping dependence of 
the chemical potential (see table (Q)). Despite a smaller 
system size we believe, that near or at zero doping, N = 8 
system shows less finite-size effects then the N — 10 when 
calculating quasiparticle properties of the system. The 
reason is that the N = 8 noninteracting fermion sys- 
tem has a six-fold degenerate level at zero energy which 
overlaps with the value of the chemical potential at zero 
doping. In contrast, the N = 10 noninteracting system 
has a large gap at [i. The scaling is therefore less ob- 
vious for the N = 10 system size, however locations of 
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the peaks nevertheless approximately scale with J and 
peak positions T qp approximately match the quasiparti- 
cle gaps obtained from the chemical potential curves for 
J = 1.6,2.0, and 2.5. 



3. Finite doping 

Spin and charge susceptibilities, presented in Fig. (|4|a 
at small doping S = 0.1 show similar high-temperature 
behavior as in the zero doping case. At high-temperature, 
Xs follows Curie-like 1/T behavior, predicted by the high- 
temperature expansion given by Eq. (||). As in the zero- 
doping case, susceptibility curves calculated for differ- 
ent J/t show scaling above T/J > T c /J ~ 0.6. With 
decreasing temperature Xs reaches a peak at T = T s 
where T s is close to its 8 = value. Even though the 
spin gap disappears at finite doping and zero tempera- 
ture, at finite temperature remains of the spin gap can 
be observed even at finite doping. Similar results were 
observed in one-dimensional calculations 113). With fur- 
ther decreasing of the temperature Xs first decreases and 
then sharply increases at even lower temperature. In this 
region the susceptibility curve can be fitted to a simple 
form Xs = 5C/T. This Curie- like form suggests that fi- 
nite doping produces nearly free localized spins. In con- 
trast to one-dimensional results |L4j|, we found that for 
intermediate values of the Kondo coupling J/t = 1.6, 2.0 
and 2.5 the local moment is reduced and equals C ~ 0.18. 
In the strong coupling limit, J/t = 10, the local moment 
reaches its maximum value C = 0.25. In the extreme 
low-temperature limit spin susceptibility should saturate 
either due to RKKY interaction between localized spin 
for small values of J/t or due to Kondo screening effects 
for larger J/t. This effect can be seen in the case of larger 
doping, 6 — 0.2, where for intermediate Kondo coupling, 
i.e. J/t = 1.6,2.0 and 2.5 susceptibility curves show less 
divergent behavior as in the 8 = 0.1 case. 
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FIG. 4. Spin a) and charge b) susceptibilities JXs, JXc 
vs. T/J at 8 — 0.1 and 0.2. Legends, given in b) apply also 
for a) and both insets. 

The charge susceptibility at S = 0.1, presented in 
Fig. (|]b), follows 1/T behavior at high-temperature, 
predicted by the high-temperature expansion given by 
Eq. (j^). Susceptibility reaches a local maximum around 
T = Tqp. For small doping T qp overlaps with its value at 
5 = 0. At even smaller temperature we observe a sharp 
increase in Xc- At zero-temperature charge susceptibility 
should equal the density of states at the Fermi energy 
which is further proportional to the quasiparticle mass. 
The sharp increase of Xc therefore suggests that quasi- 
particles are massive. 



C. Specific heat 

Results for the specific heat c v are shown in Fig. (||) 
for various values of J/t. In Fig. (|^a) we show c v as a 
function of T/t. At J/t = 0, c v has a peak at finite 
temperatures which originates from the specific heat of 
free conduction electrons. The contribution of localized 
noninteracting spins is nonanalytic and proportional to 
T8(T) where S(x) is a delta- function. At small values of 
Kondo coupling, e.g. J/t = 0.8, 1.2, 1.6 and 2.0 we ob- 
serve two peaks in the specific heat. The low-temperature 
peak that has emerged from the nonanalytic function at 
J/t = is a contribution of the spin excitations. With 
increasing J/t it shifts towards higher temperatures and 
broadens. The broadening is due to the spin excitation 
spectrum that has a bandwidth of some effective J e ff 
which increases with the strength of the Kondo coupling. 
The broad peak, that originated from the free electron 
band at J/t = 0, shifts towards higher temperatures 
with increasing J/t and becomes even broader. This is 
due to the interplay of band effects and the charge gap 
that develops with increasing J/t. The two peaks above 
J/t = 2.5 merge into a single peak which in the strong 
coupling limit scales linearly with J. At J/t = 10 spe- 
cific heat closely follows analytical prediction calculated 
within the atomic limit as seen in Fig. (^b) where c v is 
presented as a function of T/J. Our results for the two- 
dimensional lattice are in many respects similar to results 
for the one-dimensional lattice obtained by the DMRG 
method OJTl. 
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FIG. 5. Specific heat at 8 — c„ a) vs. T/t and b) vs. 
T/J calculated on a system of iV = 10. The faint dashed line 
represents the free electron result a) and result obtained in 
the atomic limit b). 

At finite doping, 8 = 0.1 shown in Fig. (|]a), small 
peaks due to the spin gap are still visible at small temper- 
atures. At even larger doping, 8 — 0.2 shown in Fig. (|j]b), 
new peaks emerge at very small temperatures (see also 
the inset of Fig. (||b)). This peaks shift towards larger 
temperatures with increasing Kondo coupling. The shift 
is approximately quadratic in J. We believe that the 
same low energy scale that gives rise to these peaks is 
responsible for saturation of spin susceptibility seen in 
the intermediate coupling J ft and 8 = 0.2. 
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FIG. 6. Specific heat c v T/t at a) 8 = 0.1 and b) 8 = 0.2 
calculated on a system of N = 10. The inset represents the 
expanded low-temperature region in b). 



III. CONCLUSIONS 

We have investigated finite-temperature properties of 
the Kondo lattice model on small square lattices. The 
chemical potential at low temperatures shows nonana- 
lytic behavior as a function of doping. The jump in the 
chemical potential is a consequence of the quasiparticle 
gap at zero doping. The quasiparticle gap scales approx- 
imately linearly with the Kondo coupling, A qp ~ 0.3 J in 
the intermediate coupling regime, i.e. 1.6 < J/t < 2.5. 
Similar scaling was recently found by quantum Monte 
Carlo simulations |L9| which furthermore show that scal- 
ing persists even below the transition to the AFM state, 
i.e. J/t ~ 1.4. In the strong coupling regime the quasi- 
particle gap again scales linearly with J as A qp ~ 0.75J 
which is not too surprising since in this case physics be- 
comes local and J is then the only energy scale in the 
system. The crossover between the two linear regimes 
occurs when the spin gap overcomes the quasiparticle 
gap, . i.e. at J/t > 2.5. Interestingly, calculations in 
one dimension show two distinct linear regimes for the 
charge gap vs. J/t fiofl . 

Two temperature scales govern the temperature de- 
pendence of spin susceptibility in the intermediate cou- 
pling regime. One scale is given by T c /J ~ 0.6 above 
which we find almost perfect scaling of curves calculated 
for a wide range of coupling strengths J. One possible 
explanation of this unusual scaling is that T c is governed 
by the charge gap A c . We have shown in the previous 
section that the quasiparticle gap at intermediate cou- 
pling scales linearly with J as A qp ~ 0.3J. Assuming 
A c = 2A qp we get a linear scaling for the charge gap, 
A c ~ 0.6J, which agrees with the value for T c . At lower 
temperatures, T ~ T g , physics at zero doping is gov- 
erned by the spin gap A s . This remains true as long as 
A s < A qp . Near the strong coupling regime the opposite 
becomes true, then the quasiparticle gap determines the 
low-temperature physics of the spin susceptibility. 

The charge susceptibility shows substantially different 
behavior than the spin susceptibility. On a smaller sys- 
tem of TV = 8 we found almost perfect scaling with J 
in the whole temperature range within the intermediate 
Kondo coupling region. This result suggests that a single 
energy scale, identified as a quasiparticle gap A qp ~ 0.3 J, 
governs the physics of the charge response of the system. 

At finite doping we found a new energy scale. It is 
reflected in the saturation of divergent 1 /T behavior in 
Xs{T) at small temperatures and in the appearance of 
low-temperature peaks in the specific heat Cy at 5 = 0.2. 
Approximate quadratic scaling of the position of the 
peaks in c v with the Kondo coupling strength suggests, 
that this energy scale could be attributed to the RKKY 
interaction between uncompensated spins at finite dop- 
ing. 
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TABLE I. The quasiparticle gap Aqp/J and- peak posi- 
tions Tqp/J and T s /J both presented in units of the Kondo 
coupling J. Quasiparticle gaps were obtained from the limit 
fi(S — > 0, T — » 0). Spin susceptibilities Xs can be in the tem- 
perature interval T/J < 0.6 and for J/t < 2.5 well fitted 
to a simple form Xs = C/Texp( — T a /T) where the effective 
moment C varies from 0.25 for small J/t to 0.5 at large J/t. 
Estimated errors where not otherwise specified are within 5%. 



J/t (N = 10) 


A qp /J 


Tqp/J 


Ts/J 


1.2 


0.40 


0.49 


0.09 ± 0.04 


1.6 


0.30 


0.35 


0.15 ± 0.03 


2.0 


0.29 


0.32 


0.23 ± 0.03 


2.5 


0.31 


0.32 


0.30 ± 0.03 


10.0 


0.56 


0.38 


0.47 ± 0.04 



J/t (N = 8) 


A qp /J 


Tqp/J 


T s /J 


[6] 


1.2 


0.27 


0.31 


0.08 ± 0.04 




1.6 


0.28 


0.30 


0.13 ± 0.03 


[7] 


2.0 


0.30 


0.31 


0.18 ± 0.03 




2.5 


0.32 


0.31 


0.31 ± 0.03 


[8] 
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0.56 


0.38 


0.47 ± 0.04 
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